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Abstract 

The aim of this article is to build on the use of tools from computational algebra initiated in 
[3], for the study of general kinetic systems, which have a wide rangc of applications in chemistry 
and biology [T31 [Ml HSl [Hl [22] ■ We clarify the relation between the algebraic conditions that 
must be satished by the reaction constants in general (mass action) kinetics systems for the 
existence of detailed or complex balancing equilibria. The main properties of these systems have 
been set by Horn, Jackson and Eeinberg [7| |9j [101 113 [Hl [Hl HO] ■ We expect to extend our point 
of view to the study of qualitative features of the dynamical behaviour of chemical intcractions 
in molecular systems biology. 

Keywords: General Kinetics, Mass Action, Complex Balancing, Detailed Balancing. 

1 Introduction 

In this article we clarify the relation between the algebraic conditions that must be satisiied by the 
reaction constants in general mass action kinetics systems, whose main properties have been set 
by Horn, Jackson and Peinberg [71 [H [lOl [I7l [181 [13 EO], for the existence of detailed or complex 
balanced equilibria. These systems have remarkable dynamic properties and have a wide range of 
applications in chemistry and biology [13l [Ml ESl [El [22] . 

We show that a reversible Horn-Jackson generalized mass action kinetics system satislying 
Feinberg's circuit conditions is detailed balanced if and only if it is complex balanced. In other 
words, under Jormal balancing conditions for the cycles (of the underlying undirected graph) of 
the reaction graph, both notions coincide. We formulate this property in terms of the algebraic 
equations deiining the corresponding varieties in rate constant space. 

In order to illustrate some of the delinitions and concepts along the paper, we will use a reaction 
network diagram which represents a nonsequential multisite phosphorylation with two sites. This 
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network models for example the Mek-MKP3-Erk system [22\ [2| [TT]: 



K46 




K1012 



The four phosphoforms, Soq, Sio, Soi, Su, are interconverted by the kinase, E, and the phosphatase, 
F. Assuming mass action kinetics, each reaction is annotated with the corresponding rate constants, 
indicated by a choice of numbering of the 14 complexes in the network. Although K32, 1^42, f^75, 
/^76) ^^109! «^119) '^1412, ^1413 j are generally taken to be very small and so the corresponding reactions 
are usually omitted, we will not ignore them in this example because we are interested in special 
properties of the reaction constants in reversible networks. 

In general, a chemical reaction network is a finite directed graph whose vertices are labeled by 
monomials and whose edges are labeled by parameters. The digraph is denoted G = (y,E), with 
vertex set V = 1,2,... ,n and edge set E C {{i,j) ^ V x V : i ^ j} oi cardinal e. The node i 
of G represents the i-ih chemical complex and is labeled with the monomial c^* = c^^^^c^^ ■ ■ ■ c|'°. 
The unknowns ci,C2, ■ ■ ■ ,Cs represent the concentrations of the s species in the network. The nx s- 
matrix of non-negative integers Y = {uij) contains the stoichiometric coefficients. For instance, in 
the reaction diagram ([T]), the node corresponding to the hrst complex E + Sqo is labeled by the 
product of the corresponding concentrations cecs^^^ In fact, we will name the 12 concentrations as 
ci,... ,ci2 according to the lollowing order of the species: 5oo, 5*10, 5*01, Su, E, F, ESoo, ESu, 
ESio, PSio, Eoi, ESoi^ So, ceCsq(, = c^c\. The rows of the stoichiometric matrix Y G {0,1}^^^-'^^ 
are ordered according to the numbering of the complexes which is reHected in the names of the rate 
constants. 

For general systems without a particular (bio)chemical interpretation, the name of the species 
will just be given by a numbering {1, . . . , s} and we record the monomial labels as the entries in 
the row vector 

^(c) = {cy\cy^,...,cy^). 

In our example, ^{c) = (c^ci, C7, C5C2, C5C3, cg, cn, C5C4, c^c^, cg, c^c^, c^cs, cio, C12, ceci). 

Each directed edge (i,j) € is labeled by a positive parameter Kj-,- which represents the rate 
constant in the reaction from the i-th chemical complex to the j-th chemical complex. Note that if 
there is an edge from i to j and an edge from j to i then we have two unknowns Kij and Kji. 

Such system will be denoted by G = (y, E, k, Y), with k = {Kij)(^ij-^^E the vector in M^q of the 
rate constants. By abuse of notation, we will just write a chemical reaction network as a digraph 
G = {V, E). Also, we will assume that a numbering of E is chosen and we will consider k G IR>o- 



2 



Let denote the negative of the Laplacian of G. Hence is the n x n-matrix whose off- 
diagonal entries are the Kij and whose row sums are zero. Mass action kinetics specitied by the 
digraph G is the dynamical system 

dc 

- = ^{c)A.Y, (2) 

where c{t) = {ci{t), . . . ,Cs{t)). 

Definition 1.1. A complex balanced system is a dynamical system ^ for which the algebraic 
equations ^{c)Afi = admit a strictly positive solution co € Mf^Q . Such a solution co is a steady 
state of the system, i.e., the s coordinates of ^(co)^k^ vanish. 

Remark 1.1. Clearly, a system ([2]) being complex balanced depend on both the digraph G and 
the rate constants Hij. A main property of complex balanced systems is that all strictly positive 
steady states c satisfy ^(c)^^ = 0. They are quasi-thermostatic [20], which in the terminology of 
[3] means that the positive steady state variety is toric. 

From now on we will assume that the digraph G is reversible, i.e. if {i,j) S E, then {j,i) G E. 
We can thus identify G = {V,E) with the underlying undirected graph G = {V,E), where E = 
{{hj}--{i,j)^E}. 

Definition 1.2. A detailed balanced system is a dynamical system ([2]) for which the following 
algebraic equations admit a strictly positive solution co € IK>o- 

- Kijcl' + Kjicl' = 0, for all {i,j} G Ĕ. (3) 

As it is for complex balanced systems, the condition of being detailed balanced depends on the 
graph G and the constants Kij. 

Note that decomposes as the sum of n x n matrices Aj^''^ for each undirected edge {i,j} € E 
of the graph G, where in rows i,j and columns i,j the matrix Aj^'''^ equals 



f^ji f^ji 

and all other entries of the matrix a]^'''^ are 0. Since the algebraic equation —KijC^ + KjiC^^ = 

means that ^{co)Ai''^^ = 0, we see that every detailed balanced system is also complex balanced. 
The converse is not true in general. Again, a main property of a detailed balanced system is that 
all of its positive steady states c satisfy —KijC^' + KjiC^^ = for {i,j} G E. 

For instance, in our example ([l]), for any choice of k G M>o the system with the following rate 
constants is complex balanced but not detailed balanced: 

Kl2 = K24 = K46 = «53 = «67 = «89 = «910 = «1012 = «1214 = «1311 = « 

«32 = «42 = «75 = «76 = «109 = «119 = «1412 = «1413 = jK (4) 
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23 ^ ^ — 47 ^ ^ _ 69 



K. 



For any a £ M>o, co,o = «(23,17,11,47,1,2,4,8,14,11,13,16) is a positive steady state of the 
system for which ^'(co^a)^^ = 0, and hence the system is complex balanced. On the other side, for 
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this choice of rate constants the system is not detailed balanced since ^ does not hold for instance 
for i = 4, j = 6, that is, for the pair of reactions 

K 

E + 5oi ^ ESqi. 

^K 

We now recall Feinberg's circuit conditions [9j, also known as Wegscheider's condition [23]. Since 
they correspond to hnear relations which only depend on the structure of the reaction graph and 
not on the particular complexes, we prefer to name them as Jormal conditions. For every cycle C 
in G, we will choose one direction and detine as the cycle in G in that direction. C~ wih be the 
cycle in the opposite direction. Although the directions are arbitrarily chosen, we will not worry 
about that since we wiU only need to distinguish between the two of them. 

Definition 1.3. A formally balanced system is a dynamical system ([2]) for which the following 
condition holds for every cycle C of G: 

n '^ii = n '^i^- 

in C+ {j,i) in C~ 

We will talk about formally balanced systems although this definition can be applied to any 
digraph whose edges are reversible and labeled by constants Kij. 
In our example ([I]), we can consider the cycle C: 



C5C3 Cii 




C5C2 Cg 



As dS]) is not satished for C, the system is not formally balanced. 
Equations ([5]) show that the set 

FBy = {k = : G = {V, E, k, Y) is Jormally balanced} 

is an algebraic variety in MS,q. 
It follows from [9j that the set 

'DBy = {k = : G = {V, E, k, Y) is detailed balanced} 

is also an algebraic variety in M|.q (see Lemma[2]). 
In turn, it follows from [3, Section 2] that the set 

CBy = {k = {i^ij){i^j)^E '■ G = {y,E,K,Y) is complex balanced} 

is a third algebraic subvariety of M>q (see Proposition [1]) , called the moduli space of toric dynamical 
systems in [3]. 

As we have already remarked, T>By ^ CBy- In fact, the main Theorem in [0] shows that 
VBy C TBy. 

In this paper we prove the following result for a mass action kinetics dynamical system associated 
to a reversible chemical reaction network G = {V, E, k, Y): 
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Theorem 1.1. Under the assumption of Jormal balancing, a reversible system is detailed balanced 
if and only if it is complex balanced. That is, 



(7) 



Our result generalizes two particular situations in which it is known that the notions of detailed 
and complex balancing coincide: the case in which G has no cycles, and the case of dehciency zero 
networks for which "DBy = TBy ([9], see also Lemma[3]below). 

Our algebraic approach tollows the lines of [3] . Our arguments easily imply that ([7|) holds at the 
level of ideals (which are radical). We refrain from giving a more algebraic formulation since it is 
straightforward and our main concern is to clarify these notions in the framework of general mass 
action kinetics networks. In Section[2]we recall known results, mainly from [313], that we state in 
a language adapted to our setting. In particular, we introduce new quotient variables which allow 
us to organize the proof of Theorem 11.11 in Section [3] 

In Section[31 following a suggestion of Martin Feinberg, we translate Theorem ll.ll to the setting 
of general kinetic systems in Theorem l4.1l and we express another necessary and sufficient condition 
for a complex balanced system to be detailed balanced in Proposition [3] 

2 Preliminaries 

Given a chemical reaction network G = {V, E, k, Y), we will denote by G the associated undirected 
graph. As we assume that G is reversible, there is no loss of intormation in passing to G. 

2.1 The minors of 

Let G be a reversible digraph corresponding to a chemical reaction network and call Gt = (yt,Et), 
t = !,...,£, the connected components of G. Up to renumbering, we can assume = A^^G) 
is block diagonal, with diagonal blocks the corresponding matrices Ai^{Gt) for the components 
Gi, . . . ,Gi. Following [3], we introduce the following detinition: 

Delinition 2.1. Consider any directed subgraph T = {V{T), E{T)) of G with n — i edges whose 
underlying undirected graph is a spanning forest of the the underlying undirected graph of G. Fix 
a connected component Gt of G and write kJ' for the product of the i^Vt — 1 rate constants which 
correspond to all edge labels of the edges in E{T) n E^. Let i be one of the nodes of Gj. The 
directed tree obtained by the restriction Tt of T to Gt is called an i-tree if the node i is its unique 
sink, i.e., all edges are directed towards node i. We will write for the product of the #Vt — 1 
rate constants which correspond to all edge labels of the edges of T^. We introduce the tollowing 
constants, which are polynomials in the (Kij): 



Tt an i—tree 

Note that each is a nonempty sum of positive terms, because as Gt is strongly connected, there 
exists at least one z-tree for every vertex i and each Kuv > for {u,v) G Et- 

It tollows from the Matrix-Tree Theorem |26] that for any i ^ Vt, the absolute value of the 
determinant of the submatrix of A^^^Gt) obtained by deleting the i-th row and any one of the 
columns, equals Ki. This (non-zero) minor is independent (up to sign) of the choice of columns 
because the row sums of A^{Gt) are zero. Compare also with the statements in |21] . 




(8) 
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Example 2.1. We will introduce a new mathematical example only for the purpose of making 
the calculations more transparent. Let G = ({1, 2, 3, 4, 5, 6}, i?, k, y ) be the following connected 
chemical reaction network: 



«12 



«25 



-J/1 



-.S/2 



«43 



«36 



«56 



-ya 



«63 



with E = {(1, 2), (1,4), (2, 3), (2, 5), (3, 4), (3, 6), (5, 6), (2, 1), (4, 1), (3, 2), (5, 2), (4, 3), (6, 3), (6, 5)}. 
For example, Ki = Yl = '^^iKs^k-gsKhk^^ + n^iK^^KesKnK^e + /«2 1/^32^63 «^43 «^52 

T an 1—tree 

+K2lK32'^63'^43'^56 + K^lK^^Ke^K^^K^l + K21K52K65K32K43 + K21K52K65K34K41 + K^lK^^Ke^K^eK^i 
+K2lK52KQnK3QKi3 + /t^l^^^^/^^^^^^l/^S^ + K 41^63^34^21^56 + «^41 /« 63 «^34 «^23 52 + /t^l^^^^/^^^^t^^/^^G 
+ «;4l/t63K34'«25'^56 + '^^l^^S^'^^^'^^^'^^^- 



2.2 The linear relations 

We now recall some results from linear algebra, which in a different language are all contained in 
[9]. Choose a numbering of the set of reactions, that is of the set of edges E oi G and form the 
signed incidence matrix Cg £ {— IjO, 1}"^^ whose column associated to the reaction has a — 1 
on row a 1 on row j and all other entries equal to 0. We denote by 

N = kcTiiY' ■ Cg) C Z^ 

the kernel over Z of the product matrix ■ Cg- 

Clearly, the vector with a 1 on its (i,j)-th coordinate and on its (j, i)-th coordinate lies in N. 
We could instead choose one direction for each pair of reactions {i,j),{j,i) in any way to get a 
directed subgraph C of G, and consider the associated signed incidence matrix Cg', with integer 

e 

kernel A'^ C Z2, since we can clearly reconstruct A'^ from A'^ and vice versa. 

The following combinatorial arguments go back to Kirchoff. We can distinguish the following 
sublattice N[ of N' . It is the Z-module spanned by the cycles of the underlying undirected graph 
G. More precisely, given any oriented cycle C we form the vector vc G {—1,0,1} 2 whose {i,j) 
coordinate equals 1 if the edge {i,j) € G' is in C, —1 if instead the edge {j,i) lies in C, and if 
neither of the edges {i,j), {j,i) is in C. 

The rank of N[ equals | — n + and a basis is formed by the Jundamental cycles associated to a 
choice of a spanning forest T of G. The fundamental cycles associated to T are those (undirected) 
cycles which are created when we add an edge in the associated undirected graph T between any 
two vertices in the same connected component of G. Note that although the number of fundamental 
cycles in a graph is fixed, the cycles that become fundamental change with the spanning forest. 

If we fix a spanning forest T of G, we can moreover choose a direct complement N^ of N[ in N' 
as follows. Consider all vectors v = {vij, {i,j) G E{G')) in N' such that Vij 7^ ^ {i^j} ^ ^(^)- 
Call A'^^ the Z-span of all these vectors v with support contained in E{T). Then 

N' = N[®N!2. 
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The concept of deficiency (5 of a chemical reaction network has been introduced and studied by M. 
Feinberg in a series of papers [H |9]. With our notations, the deficiency of the network G equals 
5 = n — dim S — i, where S is the stoichiometric linear space defined by 

S = span{yj - yj, G E}. 

Thus, dim S = rank(y* • Cc) = rank(y* • Cg). As dim 5 = § - rank{N'), we get that rank(iV^) = 6, 
so that A'^^ = if and only ii 6 = 0, and for 5 > we could choose a system of 6 generators of N^. 

In a similar way, we can decompose N as N = Nq (B Ni (B N^, where Nq is the lattice of rank 
I spanned by the 0, 1 vectors in which express the fact that the (i, j)-th column of ■ Cg is 
minus its (j, i)-th column, and A^; for i = 1,2 is isomorphic to N[ (we simply add coordinates for 
the entries corresponding to the edges not in C). 

2.3 Complex and detailed balanced systems 

We will now present some equivalencies for complex and detailed balanced systems. We will use 
the following basic result [6j: 

Lemma 1. Let k 6e a jield and ai, . . . , am £ Given a vector z = {zi, . . . , Zm) G (k — {0})™, 
there exists x = {xi, . . . , Xs) € (k — {0})'^ such that Zi = for all i = 1, . . . ,m if and only if z^ = 1 
for all A G 1/^ such that YllLi ^iO-i = 0. 

Here, we use the standard notation = Jli^i ^i' • When k = M, an easy proof can be given by 
taking logarithms. 

We now introduce new variables which are suitable for our formulation. 

Definition 2.2. Let G = {V, E, k, Y) be a reversible chemical reaction network deiining a dynamical 
system as in ([2]). For each {i,j) & E we dehne 

Qij = — Qij = -jT-- (9) 

Kji -ii-i 

Remark 2.1. The following equations hold 

qijqji = QijQji = l, forall {i,j) G E. 
We define Qij by the same formula for any pair i,j in 1, . . . ,n and then 

QijQjk = Qik, for ah /c e {1,... ,n}. 

We easily have 

Lemma 2. A chemical reaction network, G = iy, E, k, Y) is detailed balanced if and only if 

q^ = 1, for aU X e N. (10) 
Here q denotes the vector q = {qij)(ij)<^E- 

Proof. Clearly, a positive vector ci satishes a binomial equation —KijC^ + KjiC^ = if and only if 
Cq' = (lij- The result follows from Lemma[T]for m = e and {ai, . . . , am} = {Ui — yj, {hj) £ E}. □ 

We can also translate in a straightforward tashion our dehnition of tormal balancing. 
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Lemma 3. Given a chemical reaction network, G = {y^E^K,Y), the foUowing statements are 
equivalent: 

(i) The associated system is Jormally balanced, 

(ii) For every cycle C of G it holds that 

n = 1' (11) 

{i,j) in C+ 

(iii) The vector q = {Qij)(i,j)£E verifies 

q^ = 1, for all A G A^i. (12) 

Then, a formally balanced system G is detailed balanced if and only if Equations (jlOp hold for 
all A in a set of generators of N^. 

Using the results in [3, Section 2] we also have: 

Proposition 1. A chemical reaction network, G = {V, E, k,Y) is complex balanced if and only if 

Q^ = 1, for all A G N^. (13) 

Here, Q denotes the vector Q = {Qij)(ij)£E- 

Proof. We hrst claim that a network G = {Y, E, k, Y) dehnes a complex balanced system if and 
only if there exists a positive vector co € such that the following binomial equations are satisfied 

K^cl' - Kjcl' = 0, torall {i,j) G E. (14) 

It is easy to see, as in Remark 12. H that Equations ([T^ are equivalent to 

Kicl' = Kjcl' , for all {i, j) G Gt, for some t = !,...,£. (15) 

To prove that complex balancing is equivalent to (jlSp . we form as in [3j the following binomial 
ideals in Q[c] := Q[ci, . . . ,Cs]: 

I = /i + ... + /,, It = {Kicyi -Kjcy\{i,j)e Et),t = !,...,£. (16) 

We moreover dehne the ideals Tg as the saturation 

Tg = {I : (ciC2 . . . c,)°°) = {p£ Q[c] :3ueZ>o such that p(ciC2 . . . c,)" € /}. 

We denote by y>o(/) the positive variety of /, that is, the zeros of / in (IR>o)*, and similarly for 
other ideals. As Tg = {I : (cic^ . . . Cs)°°) = (/i : (cic^ . . . Cs)°°) + ••• + (/£ : (cic^ . . . Cs)°°), we deduce 
from display (8) in [3j that V>o(/g) = {c € M>o : ^'(c)^^ = 0}. But a point x with all non zero 
coordinates is annihilated by Tg if and only if it is annihilated by /. We then have that 

F>o(/) = {c e M?,o : "^{c^A^ = 0}, 

and so the system G is complex balanced if and only if there exists a positive vector co satisfying 
equations (fH|) . 

Now, we argue as in the proof of LemmaO These equations are equivalent to Cq' = Qij, for 
all {i,j) G E. By Lemma[T]for m = e and {oi, . . . , am} = {yi — yj, {i,j) £ E}, these conditions are 
in turn equivalent to 

Q^ = 1, for ah XeN. 

But from the dehnition of the vector Q, it is clear that the equalities Q^ = 1 always hold for any 
A G A^o U A^i, and so the result follows. □ 
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Remark 2.2. Note that in fact it is enough to check equahties ()13p for A in a basis of N^. For 
instance, the rank of A'^^ in the network ([T|) is 3. It is straighttorward to check that for any choice 
of constants as in ([4]) : 

Ql2 ^ Q24 ^ Qlll3 ^ Ql314 = K4K14 ~ ^ 
Ql2 ^ Q23 ^ Ql012 ^ Ql214 = = 1 

^35 ^ ^57 ^ ^89 ^ ^910 = ^t^iO ~ 

which proves again that the system is complex balancing (without needing to show a complex 
balancing steady state). 



3 Proof of Theorem 11.11 

Consider a formally balanced reversible chemical reaction network G = {V, E, k,Y). By Lemmas[2] 
and[3]and Proposition [Tl we need to show that if the constants qij satisfy Equations pT]) . then 

Q^ = 1, for aU X £ N. 

if and only if 

q^ = 1, for aU X e N. 

These relations possibly involve constants associated to edges in several connected components 
of G. In fact, it holds that modulo the tormal balancing relations, an algebraic dependency relation 
P{K) = among the (invertible) variables Qij holds for a polynomial P in e variables if and only 
if the "same" algebraic relation P{q) = is true for the (invertible) variables qij. This follows 
immediately from the following Proposition. 

Proposition 2. Let G = (y, E, k, Y) be a reversible mass action kinetics system which is /ormally 
halanced. Then, 

Qij = qij, foran {i,j) € E. (17) 



Proof. Since Equations pip relate variables quv for (n, v) in a single connected component of G, 
and given {i,j) € E then i,j belong to the same component, we can assume G is connected. 

Fix {i,j) € E. We dehne a bijection between the set of j-trees and the set of i-trees as tollows 
(see Example l3.1l for an illustration). Let T be any j-tree. 

(i) If the edge {i,j) € E{T), then let T' be the tree obtained by replacing {i,j) by the opposite 
edge {j,i). 

(ii) If the edge {i,j) ^ E{T), let Cij be the undirected fundamental cycle which is created in T by 
adding the edge {i,j). Call the corresponding oriented cycle which contains {i,j). Then, 
let T' be the tree obtained by giving to the edges of T which "lie" on Cij the direction in Cf- 
(that is, we "reverse" all these edges in T). 
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It is straightforward to check that in both cases T' is in fact an i-tree and that the map T ^T' 
is a bijection. So, we have estabhshed a bijection between the terms in Ki and the terms in Kj. 

Let T be a j-tree. We now compare the term in Kj with the corresponding term in Ki. 
If (i, j) € E(T), we clearly have that 



If instead we have that (i, j) ^ T then 



= I 

\{u,v)&C+,(u,v)i^(i,j) 

Now, by the assumption of formal balancing, we have that Y\ quv = 1 and so 

{n,v)ec+ 



n 



{u,v)£C+j,(u,v)^(i,j) 



Ki 



Therefore, 

Qij = = Qij, 
as wanted. 

Example 3.1 (Example 12.11 continued). Choose the following 1-tree T: 



□ 



K41 



«23 



«56 



K34 



«63 



Let (i, j) = (4, 1). It is clear that by reversing the edge (4, 1) G E{T) one gets a 4-tree. 
Let now (i,j) 
fundamental cycle: 



Let now (i,j) = (2,1), which does not lie in E{T), and be the corresponding oriented 





«21 




«14 


'-12 


«32 




«43 





-S/4 



-.S/3 



Then, reversing the arrows in the cycle gives the lollowing 2-tree T' 

C^l c2/2 {.2/5 



«32 



•>2/4 



-2/3 



-2/6 



«63 
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4 General kinetic systems 



In this section we generahze Theorem ll.ll to non necessarily mass-action kinetic systems in the sense 
of [8], see also [Ml Section 2]. 

The general setting is again a set of s species, with ci, c^, . . . , representing their concentrations; 
a set of n complexes, and a (numbered) set of reactions E between these complexes. We also consider 
a finite directed graph G = G(y, E, Y) whose vertices are labeled by the complexes and we keep 
our former notations. 

We assume there is a non-negative continuous real-valued rate function TZij (c) for each reaction 
(z, j) in the network, with the property that TZij{c) = if and only if c^ = for any k in the support 
of Ui. We record this information in the notation as G = G{V, E,TZ,Y). In a mass-action kinetics 
system, we simply have TZij{c) = KjjC^', with Kij £ M>o, and in this case the notations G{V, E, k, Y) 
and G{V, E,TZ,Y) refer to the same system. 

The instantaneous rate of change of Cj is given by: 

Yl yji^ji^^) - ^ yiiTZii{c). (18) 

The generalized system of differential equations that describe these dynamics could be written as 

'i=T^ChY, (19) 

where TZ is the 1 x e matrix with entries TZij, and Cq is the (transpose of) the corresponding signed 
incidence matrix we considered in subsection 12.21 

Remark 4.1. It might be useful to compare our notation with the notation in [12^ [1]. First, we 
consider a row velocity vector while it is standard to consider the transposed column vector. 
So, 

dc\ * 

— j =r7^(c)*, wherer = y*.CG G^"'''. 

Assume we have a mass-action kinetics system. We denote by K the n x e real matrix with entry 
in row indicated by complex i and column indicated by the reaction edge {i,j) equal to Kij, and 
equal to zero elsewhere. Then, TZ{c) = "^{c)K, the Laplace matrix equals = KCq and we have 

dc 

- = ^{c)A^Y = ^{c){KCh)Y = n{c)T\ 

In the notation of Y^ is called Yg, and the incidence matrices are denoted by Cg = la, K^ = Ik- 
In this general context we adapt the previous detinitions. 

Delinition 4.1. A complex balanced kinetic system is a dynamical system (jl9p associated with the 
data G{V, E,TZ,Y), for which the equations TZCq = admit a strictly positive solution co G M>o • 
Such a solution co is a steady state of the system, i.e., the s coordinates of TZCqY vanish. We call 
co a complex balancing equilibrium. 

We will again assume that the digraph G is reversible, and thus identify G = iy, E) with the 
underlying undirected graph G = {V,E), where E = {{i,j} : (i,j) € E}. 



Delinition 4.2. A detailed balanced kinetic system is a dynamical system (jl9p asociated with the 
data G{V, E,TZ,Y) for which the equations TZij{c) — TZji{c) = 0, for all {i,j} G E, admit a strictly 
positive steady state co € IK>o- We call co a detailed balancing equilibrium. 
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Again, every detailed balanced kinetic system is also complex balanced. To define formal bal- 
ancing, we need to start from a particular positive steady state: 

Definition 4.3. Given a complex balanced system at the positive steady state Cg G MS^q corre- 
sponding to the data G{V, E,Tl,Y), we say the kinetic system is /ormally balanced at cq (or that 
co is a formally balancing equilibrium) if the following condition holds for every cycle C of G: 

n ^ijico)= n ^j^(co). (20) 

in C+ in C~ 

We can now reformulate Theorem ll.il 

Tiieorem 4.1. Consider a kinetic system ()19p asociated to the data C{V, E,TZ,Y) with a complex 
balancing positive steady state cq € lR>o- have that cq is a detailed balancing eguilibrium if and 
only if it the system is formally balanced at Cq. 

Proof. Given the complex balancing steady state co G I^>o; we define constants Kij = 7^jj(co)cQ^' 
for each {i,j) € E and we consider the mass-action kinetics dynamical system ^ = ^{c)Ai^Y 
associated with G{y,E,K,Y). As 7^jj(co) = i^ijC^Q , we have ^'(co)^^ = 0, and so this new mass- 
action kinetics system is complex balanced in the previous sense. 

Moreover, as the kinetic system is formally balanced at co, we have that 

n n ^»i(co) = C W Tiji^co) = Kji, 

{i,j) in C+ {i,j) in C+ {j,i) in C~ {j,i) in C~ 

where C = Cq "-^^i^) q_ xhen, the mass-action kinetics system associated with G{V, E, k, Y) 
is formally balanced. By Theorem 11.11 it is detailed balanced. This means that every binomial 
c^j vanishes at co, implying TZij{c) — lZji{cQ) = 0, and so the kinetic system associated 
with G{V, E,TZ,Y) is detailed balanced at co. The other implication is clear. □ 

We end the paper by showing another necessary and sulhcient condition for a complex balanced 
kinetic system to be detailed balanced. 

Proposition 3. Civen a kinetic system (jl9p asociated to the data G{V, E, 1Z, Y) with a complex 
balancing positive steady state cq G K>0' following statements are equivalent: 

(i) The eguilibrium Cq is detailed balancing. 

(ii) For every cycle C in C there exists an edge i^j^j^ G ^{C) such that 7^jj(co) — TZji{cQ) = 0. 

(iii) Property (ii) holds for every basic cycle associated to any spanning forest ofC. 

Proof. The equivalence between (ii) and (iii) is clear, as well as the implication from (i) to (ii). To 
see that (iii) implies (i), let C be the digraph obtained from C by "deleting" all edges {i^, j^), j^,i^ 
in E{C), together with their corresponding labels, for all basic cycles C. Then, the associated undi- 
rected graph C has no cycles and so any positive complex balancing equilibrium co for C' is automat- 
ically also detailed balancing. Call ^^(co) (resp. ^^(co)) the Laplace matrices of the mass-action ki- 
netics system associated with G (resp. C) with reaction constants Kij = 7^jj(co)cQ for each {i,j) € 
E (resp. Kij = TZij{cQ)cQ^' for each {i,j) € E — {{i^,j^), {j^,i^), C a basic cycle of G}). But 
if co satishes the conditions in (iii), it follows that 

^'(co^A^^^co) = ^(co).4,(co) = 0. 

Therefore, co is detailed balancing for C , which together with the equalities in (iii) implies that co 
is detailed balancing for G, as wanted. □ 
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5 Conclusions 



We studied the conditions in parameter space which ensure the existence of particulary well behaved 
dynamics in general (mass-action) kinetics chemical reaction networks and we have clarihed from 
an algebraic perspective important classical notions. We plan to further apply this point of view to 
the study of biologicaUy meaningful biochemical reaction networks, in particular those associated to 
enzymatic reactions as in |15l [22l [11 [5] , where our tools from ehmination theory in the tramework of 
algebraic varieties (and in particular, toric varieties), together from results in algebraic combinatorics 
as the Matrix-Tree Theorem allow to retormulate and generaHze current approaches. 

Acknowledgments: We thank the Statistical and Apphed Mathematical Sciences Institute 
(SAMSI), USA, where this work was initiated and completed. We are very grateful to Martin 
Feinberg for his thoughtful comments. 
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